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It is proposed that there exists a class of pulsars, called weak pulsars, for which the large-scale 
magnetosphere, and hence the gamma-ray emission, are independent of the detailed pattern of 
plasma production. The weak pulsar magnetosphere and its gamma-ray emission are uniquely 
determined by just three parameters: spin, dipole, and the spin-dipole angle. We calculate this 
supposedly unique pulsar magnetosphere in the axisymmetric case. The magnetosphere is found 
to be very close to (although interestingly not fully identical with) the magnetosphere we have 
previously calculated, explaining the phenomenological success of the old calculation. 

We offer only a highly tentative proof of this “Pulsar No-Hair Theorem”. Our analytics, while 
convincing in its non-triviality, is incomplete, and counts only as a plausibility argument. Our 
numerics, while complete, is dubious. 

The plasma flow in the weak pulsar magnetosphere turns out to be even more intricate than what 
we have previously proposed: some particles, after being created near the star, move beyond the 
light cylinder and then return to the star. 


I. INTRODUCTION 

One would justifiably think that the pulsar magne¬ 
tosphere and its emission are calculable 1 only together 
with the plasma production. This is true for strong 
and dying pulsars, but (we propose) there exists a class 
of pulsars, called weak pulsars, whose magnetosphere 
and gamma-ray emission are independent of the precise 
plasma production mechanism. Weak pulsars form a 
three-parameter family, the parameters being spin (fI), 
magnetic dipole (/ jl ), and the spin-dipole angle (0). For 
simplicity we treat only the axisymmetric case, 0 = 0. 
We work in pulsar units, 

c = O = fi = 1, (1) 

and use cylindrical coordinates ( r,<p,z). 


II. WEAK PULSAR 

The pulsar is weak if it is not strong and not dying. 
The pulsar is strong if pair production rate near the light 
cylinder (the surface r = 1) is substantial, meaning > 1. 
The pulsar is dying if the plasma production rate near 
the star is not high enough, so that the proper electric 
field near the star is not fully screened. 


III. NUMERICS 

The calculated magnetosphere, Fig.l, is shown using 
the standard parametrization of the stationary axisym- 



FIG. 1: Very thick green: the boundary between the force- 
free zone (where Eo = 0) and the radiation zone. Up¬ 
per half-plane: thin black: A, integer multiples of 0.1A m i n , 
A m i n = —0.36; thicker blue and thick yellow: <p and ip, inte¬ 
ger multiples of O.lt/Wx, V’max = 4.18. Lower half-plane: thin 
blue: (r 2 + z 2 )p-, thicker red: (r 2 + z 2 )p+, isolines 0.003 x 2 k , 

k = 1,2,... 


metric electromagnetic field: 


1 The most important pulsar result of Fermi pj is the estimate 
of the median pulsar efficiency - about 15%. This shows that 
the magnetosphere cannot be calculated without calculating the 
emission, except maybe in a few cases. 


E = {—d r (p, 0, -d z cp), 

(2) 

B = - (~d z ip,A,d r ip). 

(3) 
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The calculation uses the full set of Aristotelian Electro¬ 


dynamics (AE) equations: 

B = -V x E, (4) 

E = V xB-j, (5) 

j =P+ V + (6) 

p± + v • 0±v ± ) = r, (7) 


E x B ± (B 0 B + E 0 E) 

v± =- WTT$ -• (8) 

Bl - El = B 2 - E 2 , B 0 E 0 = B ■ E, E 0 > 0. (9) 

In words: we solve Maxwell equations (4,5), with 
the current (6) from positrons and electrons of charge- 
normalized densities p± moving with velocities v±, and 
created at a rate T, (7). The velocities are given by the 
basic AE equation (8). Eq. (9) defines the proper electric 
field scalar and the proper magnetic field pseudoscalar. 
The basic AE equation and the ganrma-ray emission ac¬ 
companying this radiation-overdamped motion are fully 
derived in [2j (a) (page 2, first column, line 10). To the 
best of my knowledge, AE was first used by r 3ij. 

All details, both physics (the plasma production rate 
T, the current j inside the star, regularizations, ...) and 
numerics (grid, resolution, interpolations, ...) are given 
in the Appendix. 
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FIG. 2: Same as fig.l., with A m i n = —0.33, -i/wx = 4.19. In 
the lower half-plane: (r 2 + z 2 )p, thin blue where negative, 
thicker red where positive, isolines ±0.003 x 2 k , k = 1,2,... 

A very important result is that the full-AE magne¬ 
tosphere (Fig.l) is quantitatively close to the magneto¬ 
sphere obtained by a simplified AE calculation (Fig.2) 
which treats the plasma flow only implicitly, by adopting 


the following Ohm’s law 

. pE x B± \p\(B 0 B + E 0 E) _ 

J = - WTWi -’ P = V ' E (10) 

The agreement between the two simulations is “very im¬ 
portant” for two reasons. First, we already know [2(] (b) 
that the Ohm’s-law calculation explains the weak pulsar 
phenomenology without a single adjustable parameter. 
Second, the fact that the two methods give very simi¬ 
lar results is a strong argument in favor of the proposed 
uniqueness of the solution. 
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FIG. 3: Same as fig.l., with A m i n = —0.30, i/Wx = 4.17. 
Very thick magenta is the neutral surface - the net charge 
is positive in the equator-containing domain bounded by the 
neutral surface and the stellar surface. 

The reason for agreement between Fig.l and Fig.2 
must be as follows. From (8) we know that the true 
Ohm’s law is 

. pExB + P(B 0 B + £ 0 E) „ , 

1 = - &TeI -' p = »+ + i’- < n ) 

This coincides with the simplified Ohm’s law (10) iff the 
plasma multiplicity is zero. As we see from Fig.l, the 
plasma multiplicity is zero only in the Radiation Zone, 
meaning that the simplified Ohm’s, while exact in the 
Radiation Zone, is wrong in the Force-Free Zone. The 
agreement between Figg. 1 and 2 then means that the 
precise Ohm’s law used in the Force-Free Zone is irrel¬ 
evant. The very requirement that the Force-Free Zone 
stays force-free fixes the current in the Force-Free Zone. 

Interestingly, and also important for the analytic ar¬ 
guments of the next section, when repeated in a larger 
box, the full-AE (Fig.3) and the Ohm’s law calculations 
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FIG. 4: Same as fig. 2., with A min = —0.26, i/’max = 4.19. 

(Fig.4) do show slight differences. Now the plasma mul¬ 
tiplicity in the Radiation Zone is non-zero, and the sim¬ 
ple Ohm’s law description (10) formally fails. However, 
as we see from Fig.4, the Ohm’s law description remains 
very close to the full-AE result (Fig.3), and thus the phe¬ 
nomenological success of the Ohm’s law calculation is not 
accidental. We will say more about the plasma flow in 
the next section. 


IV. ANALYTICS 

The most striking feature of the pulsar magnetosphere 
seen in Fig.3 is the co-existence of two well-defined do¬ 
mains: (i) the Force-Free Zone, where the proper electric 
field Eq vanishes and (ii) the Radiation Zone, where the 
electromagnetic field has generic geometry with non-zero 
Eq. 

The very existence of the Force-Free Zone immediately 
raises two questions. First, for Eq = 0, Eq.(9) leaves the 
sign of Bo indefinite, meaning that one does not know 
how to calculate the velocities v±. Also, for Eo = 0, 
the charges are no longer radiation-overdamped, and one 
cannot use AE at all. 

But our numerical experiments show that once the first 
problem is fixed by some regularization of formulas (9), 
the results become independent of the regularization. We 
propose that the requirement that there be a Force-Free 
Zone bounded by the Radiation Zone (where AE does 
apply) fully fixes the magnetosphere and the densities p± 
in the Radiation Zone, while p± in the Force-Free Zone 
might depend on the regularization of (9) and on the 
plasma production rate T. Before we give our plausibility 


argument in favor of this “Pulsar No-Hair Theorem”, we 
must note another feature seen in Fig.3. 

Three surfaces, the Force-Free Zone boundary, the 
light cylinder, and the Null Surface (p_|_ = p_) intersect 
at the critical line r = 1, \z\ « 1. Inside the Force- 
Free Zone, and just outside it (except at the equator), 
the charges move along magnetic surfaces ip = const. 
If the magnetic surface intersects the Force-Free Zone 
boundary below the critical line (at low altitudes), only 
positrons flow from the Force-Free Zone into the Radia¬ 
tion Zone. If the magnetic surface intersects the Force- 
Free Zone boundary above the critical line (at high alti¬ 
tudes), both positrons and electrons flow from the Force- 
Free Zone into the Radiation Zone. 

This change in the character of the inflow from the 
Force-Free Zone into the Radiation Zone is expected. In¬ 
deed, according to Eq.(8), just outside the Force-Free 
Zone, where Eo is small, the poloidal velocities are 
v± oc B t E p ± BqB p , where B t is the toroidal magnetic 
field, and E p , B p are the poloidal electric and magnetic 
fields, Bq = B? + B p — E p . Then the product of veloci¬ 
ties is v + v- oc Ep — Bp. In the Force-Free Zone, and just 
outside it, E p = rB p (see below), giving v + V- oc r 2 — 1. 
Thus, at low altitudes, at r < 1, only one species can 
flow into the Radiation Zone, while at high altitudes, at 
r > 1, both species can flow into the Radiation Zone. 

After these observations, we can attempt a proof that 
there exists a unique AE flow in the Radiation Zone 
bounded by an arbitrary (non-AE) Force-Free Zone. The 
Force-Free Zone is described by the pulsar version of the 
Grad-Shafranov equation (Scharlemant-Wagoner equa¬ 
tion 0): 

(1 -r 2 )VV--^=-A^, (12) 

r aip 

with (f> = ip and A = A{ip). Let ip c be the magnetic 
surface containing the critical line, here defined as the 
intersection of the light cylinder and the boundary of the 
Force-Free Zone. Inside the critical magnetic surface, for 
ip > ip c , the Scharlemant-Wagoner equation is elliptic, 
and thus one can find a solution parametrized by three 
arbitrary functions: 

• the shape of the Force-Free Zone boundary 

• the current A{ip) 

• the value of ip on the Force-Free Zone boundary 

Now, ignoring for a while the magnetic surfaces ip < 
ipc we need to calculate the corresponding AE magne¬ 
tosphere outside the Force-Free Zone. It is clear that 
the result must be unique. Indeed, for the AE calcula¬ 
tion, we know the shape of the boundary, the boundary 
ip, meaning the normal magnetic field, the boundary <p, 
meaning the tangential electric field, and, as only the 
positrons may enter, we know how many positrons are 
injected across the boundary from the known A(ip). 
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Suppose the AE calculation saturates, giving some ip, 
</>, and A in the Radiation Zone. Also assume, in agree¬ 
ment with all our figures, that the isolines of these three 
fields cross the boundary of the Force-Free Zone smoothly 
(have continuous derivatives; in other words, the non- 
corotating part of the Force-Free Zone boundary is not a 
singular charge/current sheet). Then we have two match¬ 
ing conditions at the Force-Free Zone boundary. The 
smoothness of the A field follows from the smoothness 
of the ip and <p fields, because the poloidal current flows 
along the magnetic surfaces both in the force-free field 
and in the AE field with infinitesimal Eq. We therefore 
have to satisfy just two boundary conditions by choosing 
the three above-listed arbitrary functions. This seems to 
lead to an arbitrariness of the solution. We propose, in 
agreement with Fig.3, that in the elliptic region, ip > ip c , 
the position of the Force-Free Zone boundary is not arbi¬ 
trary - it must coincide with the Null Surface. Then the 
number of adjustables matches the number of the bound¬ 
ary conditions, which might lead to a unique solution. 

It would seem, however, that at high altitudes, for 
ip < ip Cl our counting argument breaks down. As we have 
seen, at high altitudes both electrons and positrons can 
flow across the boundary of the Force-Free Zone. Then 
the boundary condition for the AE calculation adds an 
extra arbitrary function - the multiplicity of the plasma 
injected into the Radiation Zone. The plasma multiplic¬ 
ity at injection, together with the three above-listed ar¬ 
bitrary functions, gives four adjustables, while we have 
only two boundary conditions. We will assume, in ap¬ 
proximate agreement with Fig.3, that the plasma multi¬ 
plicity in the force-free part of the region ip < ip c is in 
fact zero. Then the plasma multiplicity at injection into 
the Radiation Zone must be zero, reducing the number 
of adjustables to three - still one too many. But, as ex¬ 
plained by Contopoulos, Kazanas and Fendt Q, the na¬ 
ture of the Scharlemant-Wagoner equation (12) changes 
for ip < ip c . Eq.(12) provides its own boundary condition 
on the light cylinder, | d r ip = A^\ the requirement of 
smoothly crossing the light cylinder then fixes A(ip) for 
ip < ipc- This reduces the number of adjustables to two, 
equal to the number of the boundary conditions, which 
might lead to a unique solution. 

Alternatively, and also not in contradiction with the 
numerics, one can assume that the high altitudes, ip < 
ipc, are not truly force-free. There is a small Eq in this re¬ 
gion, explaining the observed absence of positrons. Then 
there really exists one more boundary of the Force-Free 
Zone within the light cylinder, the surface ip = ip c (not 
shown in the figures, while the shown boundary of the 
Force-Free Zone is in fact fictitious beyond the light cylin¬ 
der, and must be removed). The new part of the Force- 
Free Zone boundary coincides with the surface of zero 
poloidal current density — the surface of maximal value 
of | A |, emanating from the star and terminating at the 
light cylinder. If this alternative is correct, the bound¬ 
ary of the Force-Free Zone consists of the following three 
special surfaces: (i) the maximal |A| surface, (ii) the Null 


Surface, (iii) the boundary of the Corotation Zone. 

Needless to say, the above arguments are not a proof. 

V. CONCLUSION 

The weak pulsar magnetosphere calculated in this pa¬ 
per by a full AE simulation is very close to the phe¬ 
nomenologically successful magnetosphere calculated by 
a simplified procedure [2]. While some details may 
change, it appears that pulsar gamma-ray emission is ba¬ 
sically understood. 

Still, the following calculations seem doable and inter¬ 
esting 

• Using full AE in 3D, one can calculate energy- 
resolved lightcurves. Then one can measure all pa¬ 
rameters (distance, magnetic dipole moment, mo¬ 
ment of inertia, spin-dipole angle, observation an¬ 
gle) for all weak Fermi pulsars. 

• Dying (millisecond) pulsar must be a perfect tar¬ 
get for a PIC simulation (with modeled or exact 
pair production), because here the screening of the 
proper electric field is only partial, and may turn 
out to be numerically treatable. 

• Strong pulsar should be doable by adding relevant 
plasma production by photon-photon collisions to 
the AE simulation performed in this paper. 

• To really prove the Weak Pulsar No-Hair Theo¬ 
rem, and get the numerically-exact axisymmetric 
weak pulsar magnetosphere, one should solve (nu¬ 
merically, of course) the combined Grad-Shafranov- 
AE problem described at the end of §IV. Alter¬ 
natively, one may improve our bad numerics and 
convincingly demonstrate that the pure-AE calcu¬ 
lation does not depend on the regularization and 
the plasma production prescriptions. 

I thank the participants of the recent Pulsar workshop 
at Princeton, organized by Sasha Philippov and Tolya 
Spitkovsky. It was very encouraging to see how various 
different approaches close in on the unique pulsar solu¬ 
tion. 


Appendix A: Details of Numerics 

1. Grid and fields: square r-z grid; E r , j± r , B z are at 
the r-edges; E z , j± z , B r are at the z-edges; Ep±, 
j± ( j ) are at the vertices; B^ are at the faces. 

2. The star is at r 2 + z 2 < r 2 , with r s = 0.25. Inside 
the star 

P± = 0, (Al) 

j = er s (E + (ixr)xB) + j e ^, 


(A2) 
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where the conductivity of the star is large, cr s = 
200. The external toroidal current j e is chosen so 
as to give the magnetic dipole moment p = 1. 

3. Diffusive regularization: diffusion D\7 1 2 3 p± is added 
to eq.(7), with small diffusivity, D = 0.0006. Dou¬ 
bling and halving the diffusivity changes the results 
by about 10% . 

4. Velocities regularization. As written, eq.(9) gives 


B 0 = yjB 2 -E 2 + E 2 sign(E ■ B). 


We replace 


sign(E • B) —» 


E B 

y/(E-B ) 2 + aE 2 B 2 ' 


(A3) 


(A4) 


with small a = 0.001. Doubling and halving a does 
not noticeably change the results. 

5. Plasma production. We use constant plasma pro¬ 
duction rate T operating in a layer l.lr^ < r 2 +z 2 < 
1.5 r 2 , provided the altitude is above some cutoff, or 
if the proper electric field is large, Eq > aE 2 , with 
the same small a as above, just to avoid introduc¬ 
ing new parameters. The plasma production rate T, 
and the altitude cutoff are adjusted so as to min¬ 
imize the volume averaged Eq in the entire layer 
l.lrg < r 2 + z 2 < 1.5 r 2 . 


6. Interpolation. The non-diffusive part of the fluxes 
j±, namely p±v±, is calculated as follows. For each 
vertex, we calculate the interpolated electric and 
magnetic field, and use the regularized eq.(8) to 
calculate the velocity components. If, say, v+ r > 0, 
we add the current p+v+ r to the component j +r on 
the edge which is to the right of the vertex, etc. 

7. Initial condition: everything is zero. 

8. Boundary conditions: no boundary conditions are 
needed at the surface of the star. The outer bound¬ 
ary conditions are outgoing for the electromagnetic 
fields and absorbing for the densities. 


9. Resolution. Full in r and half in 2 part of the 
simulation boxes are shown. The small boxes are 
(200,400), the large boxes are (600,1200). 

We must call our numerics dubious for two reasons. 
This author can’t be impartial and after a year of 
trial and error one can get many things numerically. 
The results do show some dependence on the regu¬ 
larization and plasma production prescriptions. 
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